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ABSTRACT 


A nonlinear control system sometimes oscillate at a 
eeaency fiemenuomaneineceeral submultiple of the driving 
frequency in response to a sinusoidal input. Such a re- 
sponse is undesirable. An analog computer study 1s under- 
taken to investigate the occurrence of subharmonic oscillation 
in second order systems with nonlinearities Giana e@renizZedaby 


Baeuracion and hysteresis effects. 
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I. INTRODUCTION 


A. NONLINEAR SYSTEMS 

A system is said to be linear if its performance obeys 
Sieme ninciale sor Superposition". The "Principle of Super- 
eeoition'’ states that if C,(t) and C,(t) are the responses 
Gta System when it 1s separately subjected to two inputs 
Rj (t) and R5(t), respectively, then for all a and b, the 
mecmonse of the system to the input aR, (t)+bR, (t) is 
aC, (t)+bC,(t) and this must hold for all inputs. In other 
words, once the response to one type of signal is known, 
tie response to any other signal can be deduced directly. 

On the other hand, any system is said to be nonlinear 
fmeene “Principle of Superposition" does not hold. The 
term nonlinear systems applies to varied types of systems 
which have practically nothing in common with one another. 
As a consequence, there is no way of generalizing the re- 
sponses of nonlinear systems for a class of inputs to the 
heanolcemnor any Other inputs. The response to one type 
of input does not necessarily contain sufficient informa- 
tion about the response to some other input. Even two 
inputs which are of the same waveshape but different in 
amplitude can give responses that are entirely different. 
iimeomiaicesmthe study of nonlinear systems quite difficult. 
Even with this difficulty, there is little choice except 
to tackle the problem in some way because nonlinear sys- 


tems occupy a very important role in practical systems. 
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There is a great variety of possible nonlinear char- 
acteristics, and there are different ways by which non- 
linear systems are classified. One classification is of a 
mathematical kind based on the equations involved, however 
the study of nonlinear systems on a mathematical basis is 
Peaememely ditticult. In this study, types of nonlinearities 
encountered in servo systems are considered, namely satura- 
tion and hysteresis effects. The consideration of nonlin- 
ear elements in the operation of control systems have 
feeretrcal SitenitiGance because 1t takes linear theory nearer 


wo reality. 


IBS ¢ SUBHARMONIC OSCILLATIONS 

A linear system when subjected with a sinusoidal in- 
put produces a steady-state output which is also a sinusoid 
and oscillates at a frequency equal to the input frequency 
but may ditfer from the input in amplitude and phase. This 
is true for all frequencies. However, for a system with 
suitable nonlinear characteristic when Supe elted | bOea)S1n- 
Peoteaimnmmput, there is no guarantee that the steady-state 
Sutput is also a sinusoid. In the usual case, the output 
wave contains frequency components which are integral mul- 
meelcs Or the mput frequency. Also, the system can pro- 
duce a response which oscillates at a frequency which is 
an integral submultiple of the input frequency. This fre- 
quency response phenomenon is known as "subharmonic reso- 
nance". When the frequency of the output is 1/n times the 


driving frequency, where n is an integer, then the output 


ila 





Soomilatvon may be called the 1/n order subharmonic oscil- 
Paeron. sNhem the output oscillates with the input fre- 
G@ueney, the output oscillation is called the fundamental 
oscillation. 

The subharmonic phenomenon is of a resonant nature and 
Meeause Of this, the subharmonic component is often of large 
Tipetuce witch Causes the system output to bear little re- 
Soilpbance tO a Sinusoid of the input frequency. The output 
1s considerably different from that obtained when the sub- 
harmonic phenomenon does not occur. For systems in which 
Miomoutpiut 15 intended to reproduce the input as closely as 
moostole, aS in the motion of a loudspeaker cone or of a 
—wmVOomnlocon, this phenomenon Can cause a distortion that 
Sammot pe tolerated. in the field of feedback control sys- 
Boise it 1S interesting to know the conditions under which 
subharmonic oscillations may occur. Subharmonic oscilla- 
tions Can give false information which can be troublesome, 
SO it is desirable to eliminate the possibility of the 
Weeirmenmeer ot Slbharmonilc OSCillations. This is true in 
the case of recording seromechanisms which are frequently 
imeem COmpCELOdie inputs. Ihe possibility of the occur- 
rence of the subharmonic oscillations reduces the reliability 


of the systems. 
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Peeve Ud DESCRIBING FUNCTION 


The requirement that a nonlinearity be described by 
an equation in the frequency domain and that this equation 
Demecmpatible for use with the transfer function of the 
Jineéar components gave rise to approximating the effects 
of the nonlinear component by using a linear "approximate 
miaioter Lumection’ Or describing function. The approxima- 
mictinis Made by defining the describing function in térms 
Seure FOurler series for the component response to a sin- 
usoidal input. 

The Gesc¢cribine function is defined as the ratio of 
the neemicude of the fundamental term in the Fourier se- 
Mmiccrtor the output wave to the magnitude of the input 
Sinusoid at a phase angle which is the angle between the 
iMomsime waves and for all permissible amplitude and fre- 
Steneyeotethe input wave [Ref. 9]. In mathematical form, 
fieineminput be ASinat and the fundamental frequency term 


Hime hOUuriecr series for the output waveform 1s 
ECan omene (ie b ea Anya) 9, 


then the describing function for the nonlinear component 


1S 


Tee eae /(A,w) 
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Therefore, the describing function is but a mathemati- 
ave wemZzaevemmor a nonlinearity in the presence of a 
Senesoidal input. In analysis this sinusoid is physically 
identified as the system limit cycle and the describing 
function is used to predict its amplitude and frequency. 
The use of the describing function is based on the assump- 
tion that harmonics generated by the nonlinear characteris- 
EecSeare SUificientily filtered before being fed as input to 
Pemmonlinear element. 

Meme eCistoneOretne describing function technique is 
micwconsideration of a more Complex signal fed to the non- 
linear element. The first extension is the use of dual- 


frequency Signal, 


Vv, =a SOs eto et Db COS Mwt 


foeimput sienal to the nonlinearity with the restriction 
iitticesaebaeana nm must be real. The dual-input describing 
function can be defined as the ratio of the amplitude of 
the desired frequency component in the output waveform to 
the amplitude of the component of the same frequency in 
the input waveform. Therefore, in the determination of 
the dual-input describing function it is necessary to find 
the output waveform of the nonlinear element for an input 
composed of two sinusoidal signals and analyze the output 
to find the components of a particular frequency. In the 
systems considered here, these signals are composed of 


that component due to the input signal at a fixed frequency 
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and the ene nent due to the fed back signal which may be 
Samditterent frequency from the input signal. It is again 
assumed that the harmonics of the signal fed back are ade- 
Misc! vettleereas Ihe Concept of the describing function 
Momenus extended to allow for an input signal composed of 
two sinusoidal signals of different amplitudes and fre- 
aMeNG1es . 

Eom iustrate tne derivation of the dual-input des- 


PED Tio cnet ion, consider the cubic characteristic, 
ae a 
Let the input be of the form 
ian ci MeMioteteD Sin (Hct) 
where 


A and B are independent amplitudes, 
w and nw are any two values of frequencies, and 


8 1S an independent phase angle. 


EUbstrtutine the value of x in the nonlinear equation, the 


mesult 1S, 


[A sin nwttB sin (wt+6)]° 


“< 
It 


= A’sin®nwt+3A*Bsin*nwt sin(wtté) 


+3AB*sin nwt sin* (wt+6)+B*sin® (wtt6) 


die, 





Using trigonometric identities to express the above ex- 
pression as a series of sine waves of fundamental and 


harmonic frequencies results in 


a ie 


y = r?|3 sin nwt - 


Sin 3nwt 


= Sin(wtt6) - z sin{(2n+l)wt+6} 





+ z sin{(2n-1l)wt-6} | 


5 Sin nwt - 7 Same M2) wee ZO } 


+ 
oN) 
> 
oe) 
N 
| Paar —aennans | 
| 


- ; sin {(n-2)wt-260} | 


4 


+ 
oe, 
w 
el 
|W 


Sin (wtt6) - a Sin so (were) | 


Or 


3 
y = 7 Ce 2B) sin not - me Sone Oma 


3 


+ = (B2+2A2)sin(wt+e) - ~— sin 3(wt+0) 


4 


Sse 
4. 


+P 


| sin{(2n-1)wt-6} - sin{(2n+1)wt+90} | 


_ 3AB? 


7 | sin{(n+2)wt+26} + sin{(n-2)wt-26} | 
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Mamita peetlon of the above expression, the sum of the first 
ama third terms is the fundamental input function equation, 
mmm btiplied by the factors + (A? +2B?) and 7 (B2+2A7) . For 
the other terms, there are values of n which have to be 
considered as will be explained in succeeding paragraphs. 
Monsider the following cases: 

Clevo narmonteally related input sinusoids. In 
Mmeomease Weis 1ifrational, 1.€., n 1S not expressible as 
a ratio of two integers. In the problem being considered, 
n#41/3, 1, or 3. For the frequency nw/21, the dual-input 


deseribing function is: 


2 ACA2+2B2)sin nwt 


il 


N. (A,B 
al \PSSE White 


on 2 2 
N,(A,B) = $ (A?+#2B?) 


Mite orernequeney w/27, the dual-input describing function 


ies 2 


+ B(B?+2A7) sin(wt+é) 
pesan( were } 


i 


Np (A,B) 


ei 2 2 
Np (A,B) = zB ee) 


iMticmneswilcing dual-input describing functions, N, (A,B) and 
Np (A,B) are both independent of the phase angle, 9 but 
dependent on the input signal amplitudes A and B. 

The dual-input describing function is modified if n 


assumes a value such that the other terms in the series 


i 


. 


_ 





have Mhomeneauency Cmctiien or the input signals w/27 or 
Meyeen. 1.¢€., n = 1/3, 1, or 3. 
CPeameacruenically related input sinusoids. 
(a) n= 3. For this case there are two ad- 
ditional terms in the output waveform, y that is considered, 


imcemmtne fLOurth and eight terms, Such that the dual-input 


Pescribine function becomes 


2 
— C228. simy out - - Sim (swt+so} 


N,(A,B,3,8) = ie Sal Ree 


The above expression shows that varied phase angles 
Mee assoOciated With the terms, so the phasor representation 


iemuti lized, then 


e 2 e 
& A Une ed See 7 F— e7 Seb ttS8) 


ek 

hte) st 
fo 2 ee SIG 

N,(A,B,3,0) = = (A?+#2B2) - Zp e 


Similarly for 


° 2 ° _ 
~ (B2+2A2)e? (wtt+6) _ ~ a (wt Zo) 


Np (A,B,3,6) 7 
Bed (wtt8) 


Np (A,B,3,0) 


+ (B2+2A2) - 2b e558 


tiewresuiltang dual-input describing functions are complex 
and dependent on A, B, and 6. The resulting DIDFs for the 
Sinemet wor cases, n=l and n=1/3 as shown in the succeeding 


derivations are also complex and dependent on A, B, and 6. 
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Viemi—ion For this case, NA and Np combine to 
Give terms of frequency w/2n7. Proceeding in the same 


manner as for the case n=3, 


7 A(A2+2B2)e2%° + a (B2egaeyes (otro) 
N, (A,B,1,8) = -OoOoOororeeeem—>—” 
AeI wt 
2 . 7 ee 
= B oJ (wt 6) 4 —~ oJ (wtt28) 


+ : 
Aes #t 


2) 2 2 3B* 2 Are) OSGeo 
N,(A,B,1,0) = 2 (A2+2B2) + 2B (prs 2a2jeI9+ SB eg 


2 « 
a oJ 28 


+ 





> A(A2+2B2)eI°° 4 3B (prsza2y ed (Htt8) 
hECeBeine) = ~ 
ee ae 


2 : 2 8 
SA 5 oJ (wt 2) a: os oJ (wt+26) 
a Bed (wt ts) 


Boon 2 De 7 8 SAAD 2 Be: Bee) 20 
Np (A,B,1,8) She ONS the ) are + 7 (B 2 leat an 


3AB 56 
ie 


+ 





impr —w ly See lhisc case 15 identical with the 


ease for n=3.. 


WS 





24 (A24+2B2)e/ (wt/3) 5 = od (wt/3+6) 


N,(A,B,1/3,8) = : 
A ) Aed (wt/3) 


N,(A,B,1/3,8) = 2 (A2+2B2) - ake lt 
Ameer 56 j + 
7 i » 25 (@gesmne sd We 6) 
Np (A,B,1/3,8) es oa 
Bod (wt 6) 


eee) + - (B*+2A7) 


Np(A,B,1/3, ) 
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Peto to leM NOMEENEARITIES 


In most practical hecdbdaek Systems, the nonlinear 
emement 1S in the forward path. The basic system con- 
Sidered in this study is shown in Figure III-la, where N 
is the nonlinear element whose gain depends on the magni- 
midcwor the Input to this element and G(s) is the linear 


Pmeadency dependent portion. 


A. SATURATION 

The nonlinearities considered here are saturation, 
relay with hysteresis and backlash. Most physical systems 
Sontamimcertaim monlinearities. The most common nonlinear 
Paemonenon is characterized by a saturating or limiting 
Stmtetecristre. Limiting 15 almost universally present in 
control systems since most instrumented signals can take 
values only in a bounded range. Many error detectors, 
MicwmeacsmanGesolvyer Of synchro differential have a restrict- 
ed range of linearity. The saturation phenomenon is also 
present in amplifiers, whether electron tubes, transistors, 
Magnetic amplifiers or other devices are used. The input- 


Ceepirmeadracteristic 1S shown in Figure III-1b. 


Bye NON- IDEAL RELAY 
iiemotten Tonlinearity 1S the relay with hysteresis 
which is characterized by the curve shown in Figure I[II-lc. 


The hysteresis delays the switching operation until the 


ae. 


NONLINEAR ELEMENT LINEAR ELEMENT 





(a) Block diagram of basic system. 





(b) (c) (d) 


NONLINEAR LINEAR NONLINEAR 





(e) Basic system with two nonlinear elements. 


FIGURE III-1: System block diagrams. 
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Mapes epast the desired zero point. A system using this 
fen or relay must have a limit cycle in steady state. 

iiemrelay sis a device to apply the power which drives 
micmriogd.) Inta relay control system, the power amplifier 
is a relay device. The relay amplifier is desirable be- 
mause lt Gan be Simple, rugged, compact, and relatively 
mieape wille meeting high load power requirements. Ex- 
amples where relay control systems are typically applied 
are in space vehicle attitude control and aircraft and 


missile adaptive control systems. 


Cc. BACKLASH 

Another commonly encountered nonlinear phenomenon in 
iestemecndnical linkages is backlash. In control systems 
arenes present im most mechanical gears. tJhis is true 
vem £oOr mew ones and with use, the backlash increases. 
inesimple physical terms, backlash in gears is the free 
Space between adjacent teeth. It occurs whenever there | 
is inexact linkage between mechanical parts, i.e., when 
there is slack between mating gears. The backlash char- 
heen rct ers siown on Figure Iil-ld, where A represents 
Paesamount Of backlash. 

In one way or another, linear and nonlinear systems 
Woibinemeounter the backlash problem if they have a gear 
Mertiemeimethe tieclda Of Control systems, a system is in- 
cluded in the class of unstable systems if it continuously 


oscillates with constant amplitude and zero input (limit 


ZS 





eetle). Iitatersand Pastel [Ref. 9] show that for a second 
PmeorescrVOeWwiun Dackilasi, the existence of the limit cycle 
PreidcwOneeicmvdalice Of tne damping factor, c. For t > 0.29, 
Biemoystem m= stable so that for 5 < 0.29 the effect of 
backlash is to cause an unstable character. If the phase 
portrait is drawn for the system, the phase portrait will 
aon the backlash ettect near the origin. For systems 
Prentiss dcomred to have the phase trajectories termin- 
ate at the origin, the backlash effect causes serious - 
erouble. 

The three nonlinearities discussed above were simula- 
meoaeon the electronic analog computer. The operation of 
time simulators for the nonlinearities are explained in 
Popenadix Baseline Characteristics of the nonlinearities were 
Pmocuiced very nicely as Shown in Figure III-2. 

Pretec oe bi=4oeand T11-5 show the block dia- 
grams of the second order systems with the indicated non- 


linearity which were simulated on the analog computer. 
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(c) Relay with Hysteresis 


Settee i 2.ssoamulation output for nonlinearities. 
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FIGURE III-3a: Block diagram of second order system with 
saturation. 
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FIGURE III-3b: Analog simulation using digital logic for 
Sul acron ot Saturation. 
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FIGURE III-4a: Second order system with non-ideal relay. 
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FIGURE III-4b: Analog simulation of system with relay. 
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IV. SUBHARMONIC GENERATION 


Consmadey thespaste System Shown in Figure JJI-la with 
a nonlinear element in the forward path. If the nonlinear 
meeiient 15 replaced by an element of unity gain, the result- 
mroelinecar SyStem 1$ assumed to be stable. 

eoetecdmendieilaei.. the "Principle of Superposition" 
applies to linear systems. The response of the closed loop 
system to a given signal is independent of any other signal 
which may be applied at the same time. When a nonlinear 
Seement ts included in the feedback loop this is no longer 
MiPmecdase. Ine response of the nonlinear element to a given 
Signal is modified if another signal is applied at the same 
Emme. ) Ine €ftect of the other signal may be considered as 
wemvincgethe fain Of the nonlinear element for the primary 
mipieesmenain “Ihe case of interest 1s when the additional 
Signal input is also sinusoidal with a frequency that is 
aueilcesraliemultiple of the primary input frequency. In 
DeeitsmGasc, the OULpUtL Compcnent of the angular frequency, 
w may be modified in amplitude and phase by the presence of 
the additional input. This modifies the loop gain of the 
system. The phase change is caused by the intermodulation 
PuOdvcucmecne tated by the nonlinear Clement. 

Hiemoedpilityeorwene closed loop system is affected by 
varying the gain of the nonlinear element since the nonlin-- 


ear element is a part of the feedback loop. For a range of 
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BenuUsSoOlGal imputs, the response of the nonlinear element 
may be so modified that the closed loop is rendered oscil- 
Hacory. Where the system 1s such that the open-loop fre- 
Micniey response plotted on the complex plane do not 
intersect with the negative real axis, the nonlinear ele- 
ments have to introduce an adequate phase change in order 
Gere stable oOScillation can occur. This implies that the 
frequency of the oscillation must be related to the input 
mmeaueicy. What 15 Of interest here is where the fre- 
quency of the oscillation is an integral submultiple of 
Pee Orimary input frequency in which case subharmonic os- 


cillation is generated. 


31 





V. DISCUSSION 


iiieesystemsesnown in Figures [II-5, [II-4, and III-5 
with the indicated nonlinearities were simulated to in- 
vestigate the occurrence of subharmonic osciallations and 
bowobserve the effect of changing the variable parameters 
in the systems such as, the damping factor, z and forward 
linear gain, K. Mimics li omand Liil-5 and the amount 
Sumeibysteresis, ZC in Figure III-4, in the generation of 


Sulasarmonic oscillation. 


A. SYSTEM WITH SATURATION NONLINEARITY 
Figure III-3a shows a second order system with a sat- 
Meatvomenonlinearity. Ihe nonlinear characteristic is 


defined as: 


e =e ee Ca <7 
O al 1 

e = 2 e. > ? 
O al 

e = -2 e. < -2 
O al 


piceeenem cain Of the linear portion of the saturation is 


Mb eine mtranster function of the linear part is, 


K 


m 
G(s) ~ S(s + 0.4 


Considering the linear part, the closed loop transfer 


momction 1s 


B. _ Rn 
Op a at cat kK 


OZ 





ramen CNaracter1sStic equation 
s*¢ + 0.4s + K_ = 0 
m 


iie natural frequency, om TitdmencmGanpine factor, © 15 


determined from the following relationships: 


= 
i 
. 
= 


mY 
N 
no] O 
am I 
lc 
ee 


m 


meecehe Simulation of the system in the electronic analog 
computer , K mceveartcdsiromn lL to 100 giving a variation 
Mmemrronm 0.2 to 0,02. Typical results obtained by the 
fmlbation are shown in Figures C-1 through C-5. 

Ogata [Ref. 6] discussed an analysis of subharmonic 
responses in systems with a nonlinear element whose input- 
Output characteristic curves shows saturating or limiting 
Simactenisties and showed that for nonlinear characteris- 
miOmectinvesetiiat are Skew symmetric, odd-order subharmonic 
esemlvations Can occur. 

item: Zco-deerec criterion’ for predicting responses of 
Saturating feedback systems was introduced by Douce [Ref. 1]. 
In the complex plane, two lines are drawn from the point 
as -1 on each side of the negative real axis and each line 
Maxeccwateanglewon 28... This 15 illustrated in Figure V-1. 
The region bounded by the two lines represent the critical 
region for a system with saturation nonlinearity. Utilizing 


PMbommencmNyaduist diagram of the linear part of the system 


OS 





von 


Glyw) = 
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PrGUREMV- 1: 


The 28-Criterion (Douce [Ref. 1]. Shaded area 
Shows region where 1/3 subharmonic can be gen- 
CRatcaema systems Containing @ limiter. 
DOtEowehiiies Tepresent open loop frequency 


response of linear system with indicated value 


Of AK 2, 
m 


1.5) 





maovides sutticient information to predict possible occur- 
rence of subharmonic oscillations. Subharmonic resonance 
memrossiplec whem the Open loop frequency response of the 
minicar DOrt1OnM6enters the critical region. Figure We Ih 
Emewcethe Oper Loop frequency response of the linear part 
for different values of Ka: Also indicated is the damping 
factor. The figure shows that for Ke ae Ome @ Oo 44 en 
Gremopen Loop frequency response of the linear portion 

lies outside the critical region so that subharmonic OS- 
@m@eraeion cannot occur. For i 2a ero Lexan lc Ks = $0, 
Miemianear open loop frequency response is inside the crit- 
feareregion. lhe results of the analog computer study 
iemeecsacctirately with the '"28-degree criterion" as can be 
seen in Figure V-2 which shows the domain of 1/3 subhar- 
Meme Oscillation as a function of input frequency and 
fiipenowiaetor. sc. Ine Ligure also shows the effect of 
hMpiemamplatude, 1.€,, for preater input amplitude, the 
domain decreases and the system has to be more lightly 
damped to sustain subharmonic oscillation. 

Subharmonic oscillation of 1/3 order was the only one 
observed for the system with skew symmetric saturation non- 
linearity. Figure V-3 shows the domain of 1/3 subharmonic 
as a function of input amplitude and input frequency. It 
can be seen that the range of input frequency where sub- 
harmonic oscillation is sustained decreases with increasing 
input amplitude. The same observation can be seen in 


Figure V-2. It can be deduced that subharmonic oscillation 
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Saoccur only for a2 range of input amplitude. For a small 
mipute amplitude such that the error signal to the nonlinear 
Seemente 15 less than the saturation level, no subharmonic 
oscillation can possibly occur because the system will just 
be linear. Also shown in the figure is the variation of 
eiepue amplitude. For a particular damping factor, the 
Seamoee GG —- 0,02 in Figure V-4, the output amplitude varies 
memonown by the Curve irregardless of the input amplitude. 
itesouLput amplitude is dependent on the range of input 
mecauency that the subharmonic oscillation is "locked on". 
Figures V-4 and V-5 show the output amplitude variation 
Domi romeme values of ¢€. Ogata [Ref. 6] concluded that 
subharmonic oscillations do not occur as response of the 
GuenueE Of the nonlinear feedback control system when the 
Pept -inpwe amplitude ratio for a given amplitude of a 
Tiicoldaleainput 1s 2% 3 db or less. This is equivalent 
Boomer attowo: l.256 % J.41. The results of the simulation 
Peldyvedo mot agree with this criterion as can be seen in 
Figures V-4 and V-5 where the lowest output-input amplitude 
PaietocmunewilmeZoeand 1,125; respectively. The possibility 
of subharmonic oscillation is greater for greater output- 
Maput ratio. 

Subharmonic oscillation of 1/2 order was observed when 
the saturation nonlinearity is skew asymmetric. The sat- 
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Output amplitude, volts 


FIGURE V-5S: 





Inpue Frequency, rad/sec 


Output amplitude for 1/3 subharmonic due to 
AcE Omi apiiGceamplitude =916 volts. 
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The domain of 1/2 subharmonic is shown in Figure V-6. Fig- 
Mires o-1 through C-5S show the signal waveforms for 1/2 and 
i seorders subharmonic oscillations and also the waveforms 
of the error signals. Figure C-2 shows the build up of 
1/2 order subharmonic and Figure C-4 shows the decay of 
iWeseorder subharmonic to fundamental oscillation. 

leocwmloulcemandmlavesly |Ref. 12] presented a method 
eeeiisime the dual-input describing function technique to 
eeuway SuUbharmonic oscillation. Figure V-7 is the describ- 
ing function loci which shows the response of a saturation 
nonlinearity to a sinusoidal signal of amplitude B and of 
eeorenal at three times this frequency and of amplitude A. 
awe cetapibates other describing function loci for dif- 
femeie values of amplitude B. 

MecupmarnionLe oscillation of 1/n order where n is an 
miecoer will oscillate at a frequency that is 1/nth of the 
mimic reqiency.) [oO investigate this subharmonic, let the 


mput to the nonlinearity be 
MacOS iit B Cos (wtté) 


The signal A cos nwt is assumed to be known and the des- 
Swibim@ecimetion tor the Signal B cos (wttd) is superimposed 
on the open loop frequency response locus of the linear 
Part. The value of B, such that the describing function 
cuts the frequency response locus for a frequency, w 1s 

the amplitude that corresponds to an overall open loop gain 
of unity. So this gives a possible steady-state amplitude 


of the resultant subharmonic oscillation. 
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The amplitude and frequency of possible 1/3 order 
subharmonic oscillation is found as follows. Assume that 
the amplitude, A of the input signal is known. For the 
System to oscillate continuously at a frequency equal to 
wieectnerd Of the input frequency, the overall open loop 
Painemust De Unity. Then the conditions for continuous 
oscillations may be evaluated by means of the describing 
Pfmer rom technique. The input to the nonlinear element 
1s 


PCOS moOut I+ BR COSMUDt tS) 


Know the value of A, the describing function loci that 
Searespomd co this value of A are superimposed on the open 
loop eens response locus of the linear system. Refer- 
ing to Figures V-7a and V-7b this shows the describing 
mometion loci for A = 4 and A = 8, respectively. The open 
loop frequency response locus of the linear portion with 


Ceikransitenr finction of 


ee 50 
CO Suto 4) 


is superimposed. The point at which they pass through the 
point representing the frequency response at the subhar- 
monic frequency determines the value of B, which is the 
amplitude of the subharmonic oscillation. 

Figure V-8 shows the output amplitude for 1/3 sub- 
fommlon@enosciillation tor a second order system with sat- 
Metronmmonlinearity where the Saturation level is + l. 


The curves were obtained from the results of the analog 
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|= 
aN je (jw 0.4) 
Re ~14¢ ~J2 
ae 
Sie ee 


BiG a aa Yesermbang function loci for A = 4. 
Shows response of saturation nonlin- 
Cire eOmaesimusoidal signal of 
amplitude B and of signal AE three 
times this frequency and of ampli- 
Gude Ret. 12]. 
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PGURE  V- 7b: 


Describing tumetron Lecr for 

A = 8. Shows response of sat- 
uration nonlinearity to a sin- 
Wceohagiestonaleor amplitude 8B 
and of signal at three times 
this frequency and of amplitude 
Maroc moon requcncy Fesponse 
FociicuoOtultIneam system super- 
iMpesed =a ked. 12), | 
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Table V-3: Amplitude of system response and error signal 
for system with saturation level il. 
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Sompmter study, Also indicated are points calculated using 
the describing function loci of Figure V-7 and those tabu- 
Mmiced in Ref. 12. 

Pies@imeteulty with the dual-input describing function 
Hecimique 15 the derivation of the dual-input describing 
Mmimceron atself. if the nonlinear characteristic can be 
expressed as a simple power series or simple output-input 
relationship such as the cubic nonlinearity considered in 
@iapter Ii, the dual-input describing function can be de- 
menmined by simple algebraic analysis. However, in the 
Sse ol practical feedback systems where the nonlinearity 
1s rarely expressible as a simple power series, like in the 
@7se, 08 Saturation, relay, and backlash, the algebraic an- 
alysis of the nonlinear element becomes difficult. Graph- 
ical methods are possible but very laborious. Due to the 
complexity of these nonlinearities, the use of computers 
Pomedetermining the relevant dual-input describing func- 
te_oOm 1S necessary. 

Pomethiespurpoese Of this study, the derivation of the 
dat -input describing function for the nonlinearities are 
not included. The DIDF derivation can be treated as a 
separate problem and the results can then be applied to 
mvtommeor linear systems to study the effect of the 


montimearity. 


B. SYSTEM WITH NON-IDEAL RELAY 
The other second order system simulated to observe 


the occurrence of subharmonic oscillation contained a 
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eouplenvaltied nonlinearity. Figure I]I-4a shows a second 
order system with a relay with hysteresis. Referring to 
this figure, the nonlinearity can be defined by the fol- 


lowing equations. Assume the input to the nonlinearity 


1s 
e. = A sinwt 
1 
then 
e = -a Oe wit <6 
O 
ey = a pace wWe Sen tp 
where 
g = sin} c 
A 





and the plot of the inverse of this describing function 
for a = 1 (constant) and three different values of c, where 
Vepiomume salount Of hysteresis, is shown in Figure V-9. 
Superimposed on the -1/Gy Slmvestsmenc Dlot or Gl jw) of 
the linear portion of the system. The figure shows that 
Pm actrecular amount of hysteresis, the point of inter- 
section between the G(jw) curve and the relevant -1/G,) 
curve represents a limit cycle whose amplitude and fre- 
quency can be determined easily. For example, consider 

c = 0.2, the limit cycle frequency is approximately 1.8 
rad/sec with an amplitude of 0.46. The Bode plot of the 


Simple linear portion is shown in Figure V-10. 
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Describing function for relay with hysteresis. 
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The system is Figure III-4a has a sinusoidal input of 
known amplitude and frequency. The output of the relay, 
f(t) is assumed to lose its higher frequency components due 
to the low pass characteristics of G(s). Then, the input 


to the nonlinearity can be written as 
eos si Gta see OS Wit Gg) 


where A cos nwt is the primary input signal of known anm- 
meitude and frequency. The dual-input describing function 
of subharmonic oscillation developed by West, Douce and 
Iuvesily [Ref. 12} can be used if the dual-input describing 
fIMNEetLON 1S available. 

In the ‘analog computer study for the system, subhar- 
mene oscillation as high as 1/2 order and as low as 1/13 
dere were Observed. The subharmonic oscillation once ini- 
@uared 15 Stable and stays ‘‘locked on” a reasonable range 
fiomeredicney. Ihere 1s an overlapping of domains of dif- 
ferent subharmonic order, as shown in Figure V-ll. The 
subharmonic oscillation that is stable is the subharmonic 
Seaemeliiitiated first. It can.be seen from Figure V-11l 
that odd-order subharmonics have a much larger domain than 
even-order subharmonics and that the domain decreases for 
lower odd-order subharmonics. The figure also shows the 
iplrtudesor the Output which is the amplitude of the sub- 
iteanonmemosecillation. ThiS 1S the amplitude curve for all 
SupManmomi1e modes. If the output oscillates at a frequency 
mreeioerad) sec, then the output amplitude is 0.225 irre- 


gardless of the subharmonic mode. 
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hapes-5. Simulation Data for Non-Ideal Relay (continuation) 
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1/3 Subharmonic 


1/5 Subharmonic 








4 Vi e) sie oe 


OUTPUT FREQUENCY, rad/sec 


FIGURE V-ll: Domains for subharmonic output for system 
Wetiemetavowae— i Cc = 0.2. 


Sf 





Cw ie snows the domain of subharmonic oscilla- 
moms: Orea =) '!eand € = 0.1 adapted from the work of Sawaragi 
Pocmakashin|kRet. S|. The points plotted represent data that 
were: obtained from the analog computer experiment. 

Figures C-6 through C-16 show the simulation results 
for the system with non-ideal relay. In the system with 
Saturation nonlinearity, subharmonic generation was not 
feeeestarcing. The occurrence of subharmonic oscillation 
Paoecepenaent On-the initial conditions, such as, the sud- 
wemeenmange Of input frequency or input amplitude. For the 
eeeemewiclm non-ideal relay nonlinearity, subharmonic oscil- 
Pomeonemvias Seli Starting. This can be seen in Figure C-6 
which shows the system response changing from fundamental 
Paemebacron to 1/5 subharmonic oscillation. Figure C-6 
miseushows the relationship between the amplitude of the 
output when oscillating at the fundamental frequency and 
faecmmoscr iilating at the subharmonic frequency. Figure C-10 
shows the change from subharmonic to fundamental oscilla- 
tion, metro l/s Order subharmonic to fundamental. The 
amplitude of the output at the fundamental frequency is 
practically negligible compared with the output amplitude 
Mitem@eoscililating at a subharmonic frequency. Figures C-/7 
through C-16 show the changing modes of the subharmonic os-' 
cu wioneaue to chance im imput frequency. There were in- 
stances when the output oscillates at two different subharmonic 
modes at the same input frequency. Figure C-16 shows the 
output oscillating at two different subharmonic modes; name- 


ly, 1/11 order and 1/13 order at the same input frequency. 
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C. SYSTEM WITH TWO NONLINEARITIES 

The third system simulated to investigate the occur- 
fmomeesOr SUDharmonie oscillation contained saturation and 
friction controlled backlash nonlinearities as shown in 
Peete til-Sa. It is a two loop system with a tachometer 
feedback for the inner loop to provide the necessary damp- 
ing to make the system stable. 

powilesstandand procedure, the characteristic equation 
1S 

SO MNS 

Pietsecleéar that the damping factor, ¢ is dependent on or 
and K,, ess 


he 
ee ge a 


mM 


Han [Ref. 4] presented a procedure to predict the 
Srrstenee Of limit cycles in a system with two nonlinearities 
Mlcwmasetne system shown in Figure IlI-5a. For the purpose 
ond Bpeeivine EReCROGcell Geneemon —SUbmarion1c OoS¢cillation in 
the system, assume that the system oscillates with a limit 
cycle so that the input to the first nonlinearity is of the 


form 


eGCOsmHNuUteta bh cos (wt te ) 
where A cos nwt is the primary input signal of known ampli- 
tude and frequency. The linear frequency dependent ele- 


ment, G(s} is a low-pass filter which is assumed to adequately 
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Beremuare the higher frequencies so that harmonics can be 
neglected in Can It also follows that the input to the 
second nonlinearity is a sinusoid. 

. The analog computer simulation was performed for dif- 
ferent values of i and K, ana tie results are tabulated 
mimbable V-6 and V-7. Figures C-17 through C-23 show the 
Signal waveforms obtained. Figure C-17 and C-18 show the 
femetorms for 1/3 and 1/5 orders, respectively. It can 
meomseen that the signal en which 1s the input to the back- 
lash nonlinearity oscillates at the subharmonic frequency, 
Wace ti Figure C-17, Can 1s oscillating at a frequency 
erates 1/5 of the input frequency and in Figure C-18, at 
1/5 of the input frequency. The system output, 8. oscil- 


lates at the same frequency as e.. with a phase difference 


Tera 
of 180°. This shows that the subharmonic oscillation ie 
dieweo une effect of the saturation nonlinearity. The ef- 
Poet onetne backlash nonlinearity can be analyzed by using 
Picmconventronal describing function since the input to the 


backlash nonlinearity, e.. is only composed of one sinusoid- 


in 
al signal. Figure C-19 shows the change of fundamental os- 
Smbiation to 1/35 subharmonic oscillation and Figures C-20 
mide -o2ieshow the build up of subharmonic oscillations. 

It can be seen also from these figures that the subhar- 
monic oscillations are self starting as compared with the 
system with saturation nonlinearity only. 


The domains of subharmonic oscillation as a result of 


the simulation study are shown in Figures V-13, V-14, and 


61 





foro.) the domain of the 1/3 subharmonic oscillation is of 
the order of the natural frequency. For lower order sub- 
Meommensc modes, i.e., 1/4, 1/5, 1/6, 1/7, ., subharmonic 


Seeiletation occurs at a higher input frequency. 
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Table V-6: Domain of subharmonic for system with two 
nonlinearities, ¢ = 0.283, oO, = 14.14 
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VI. CONCLUSIONS 


The investigation of subharmonic oscillation in non- 
linear systems was undertaken by means of an analog com- 
puter study on three nonlinear second order feedback 
control systems. The following conclusions are summarized 
meom the result of this study. 

When the output of the system with saturation non- 
linearity is oscillating with the input frequency, sub- 
harmonic oscillation cannot occur unless some kind of sudden 
disturbance is given to the input. This sudden disturbance 
may be a sudden disturbance may be a sudden increase in 
mpueeamplaitude or frequency. This is not the case for 
the systems with hysteresis nonlinearity. In the system 
with non-ideal relay and backlash characteristics, sub- 
harmonic generation was self starting. 

iittemuecsibulity of the occurrence of subharmonic os- 
Geelaudiom 1s great when the output-input amplitude ratio 
is large. In the case of saturation nonlinearity, there 
is a minimum output-input amplitude ratio for subharmonic 
Oscillation to occur, i.e., the output-input amplitude 
Paetomise coreater oF equal to 1.125. For the system with 
mon-adeal relay, subharmonic oscillation can occur for 
output-input amplitude ratio less than unity. 

Hhe "28-deeree criterion" introduced by Douce [Ref. 1] 


Gime omaeccirate results in predicting the possibility of 
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subharmonic responses in a system with saturation nonlin- 
Pwerey. INas 15 a useful technique in the elimination of 
subharmonic oscillation for systems with a limiting non- 
immearity. 

Subharmonic oscillation once initiated is sustained 
over a reasonable range of frequency. There 1s an over- 
iippineg of these frequency ranges for different subhar- 
monic modes. In the region where two or more orders of 
subharmonic oscillation can occur, the order that occurs 
first dominates and tends to continue and 1s quite stable. 
The frequency of the subharmonic oscillation has a minimum 
depending on its order. The frequency for 1/3 order sub- 
harmonic is usually in the neighborhood of the natural 
Pmeeoaducmhey Or the system. Subharmonic oscillation is de- 
Pendent on the input signal amplitude. It can occur only 
Mecniieaecertain range Of signal amplitude. It can be con- 
cluded that subharmonic oscillation can occur only with a 
Gertaim range o£ input amplitude and frequency. 

The system with a limiting nonlinearity is dependent 
on the damping factor for subharmonic oscillation to occur. 
The system has to be very lightly damped, i.e., ¢ < 0.1 in 
order for the system response to oscillate at a subharmonic 
frequency. If the input amplitude is increased, the system 
has to be much more lightly damped for subharmonic oscilla- 
PiOmmranOcCur. 

The dual-input describing function technique developed 


by West, Douce, and Livesly [Ref. 12] appears to be a very 
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useful technique for predicting quantitatively the ampli- 
tude and frequency of subharmonic oscillation. Recommend- 
Peon ture Study 1s the derivation of relevant dual-input 
describing functions (two sinusoidal input) and ap lieereion 
Seene results to any linear system with the relevant non- 
mmMcarilty. Quantitative analysis of the occurrence of sub- 


marmonic oscillations should be undertaken. 
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(2) 


(3) 


(4) 


(5) 


APPENDIX A: PROGRAM SYMBOLS 


Servo Set Potentiometer 


(x): e = Ke. where K < 1 
1 O O al 


Summer Amplifier 


e. 
il , 
eS 5 (O54 + 10e,,) 


@?) 
‘e) 
Oo 


Say 


Inverter 


miter akor 


Stal 


OU = - 
ag fo os Ce a ee 


e. 
15 


Comparator 
OMtput Logic: 


Algebraic sum of X + Y > Q 


C 


C 0 volt (false) 
Algebraic sum of X + Y < 0 

Ge oO lt mera lisic) 

Gee somvolts Mrruc) 


Ove hts (true) 


m< Ps 
Or ‘-) 
" i 
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(6) 


(7) 


(8) 


Digital/Analog (D/A) Switch 


Mogremcon tio ll: 


Hiieeswmech 1s —ClLOSea if either 


a “6 Oe voOltcuis appl Ted On. the 
Piptewtswlet te Open; the switen 
to sOpcieieUeTOlt TS applied. 

OF Gave 

Ly zs 

= wey 

ay, 

AND Gate 

aa 

6 
L, + L 
oy. 
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APPENDIX B: ANALOG SIMULATION 


10 Nom linear ices 
(a) Saturation 
PeClGemeabamonowse tne Saturation Characteristic 
eve and Figure B-lb 1s the circuit diagram for the simu- 
Prone usiie the analog computer. Referring to Fig. B-1b, 
the simulation works as follows: 


(1) a = a 


Pacebuarersum of Cl < 0! > Cl false, Ci 
uve and also the algebraic sum cf CZ < 0 + C2 false, 
C2 true such that the output of the AND gate = 0 and 
D/Al is open. For Cl true, D/A2 is closed and oe M077 


mrcees D/AS 1S open, therefore, 


Wimevlcebrare sumo Cl > 0 > Cl true, Cl 
false and the algebraic sum of C2 < 0 + C2 false and C2 
true. The output of the AND Gate = -6 volts and D/Al 
momellosed, For Cl false, D/A2 is Gpeumand G2 etalse., 


Dyeaseaseopen, therefore, 


Algebraic sum of Cl > 0 + C1 true, Cl 


false and the algebraic sum of C2 < 0 + C2 false and C2 
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Paauiceeb-o la. Gaturation Characteristic curve. 


ei Se D/A 4 


D/A2 


D/A L 


= D/A 





FIGURE B-1lb: Analog simulation for saturation. The 


View rch nomtuomtrid of anplitfier 
and 100 is the reference voltage of the 
computer being used. 
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muucee ihe Output of the AND Gate = -6 volts and D/Al is 
euosed. Por €i false, D/A2Z iis open and C2 false, D/A3 


Momoepen, therefore, 


eeciomamemsunmeate Cl > 0 > Cl true, CL 
feese and the algebraic sum of C2 > 0 > C2 true and C2 
false. The output of the AND Gate = 0 and D/Al is open. 
Pomeclietalse, D/AZ 15 open and for C2 true, D/A3 is 


closed, therefore, 


In the description of operation above, -6 volts represents 
Ploered!) fyand 0 volts represents a logical 0. The logic 
control for the D/A switches is the output of the compara- 
for or of the AND Gate. 
(b) Relay with Hysteresis 

Seren > = Jae roeriic wenanacrenistie Curve OL a 
WOolterdcaiierelay with hysteresis and Figure B-2b is the 
ieaeranm tor analogs simulation. The Track-Track (T/T) unit 
is composed of two D/A switches with a common G terminal. 
Pielitmemsb=7651S the analog Circuit for the T/T unit. Switch 
T is closed and T is open when a logical 1 is applied to 
Biome puteor the input 1s left open: When the input 1s a 
logical 0, the switch states are reversed. When the 


Switches are open, the resistor switch ends are grounded. 
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micUee B-2a: Characteristic curve of relay with hysteresis. 


2 
aa 


G Ss; 

+B T 
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_ ! 
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i 

* ces 
c/a = 


MicUREeb-= 2b: Analog simulation of relay with hysteresis. 





mueimeme= 2c: Analog Circuit of track-track (T/T) unit. 
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fe) ee backlash 
Tiembacwlascmmeenatacter1Stic Curve is shown in 
Feune B-3a and the analog simulation circuit diagram is 
Shown in Figure B-3b. The backlash Braracrer is cac 1s 
achieved as follows: 
PoOmeeminencedces irom zero nothing happens to eg 


1 


mci 1 (e. owe wc minemic that point the integrator 


O 

output will begin to change at the rate m(e; = Cy) > and 

ey will approach (e. = a) through the feedback correction 
) 


becomes less than a resulting in e'' going to zero and the 


Gaetne loop. AS soon as eC. ehangesmailrection, (e. cao. 


output of the integrator holding its value. As e; con- 
eimues to decrease, eC remains unchanged until e' becomes 
ocomtianeby at which time, the output of the integrator 
Bemins tO Change at the rate m(e, : Sed ¢ Again when e. 
eimecs arltrection, e'' poes to zero and the output of the 
integrator remains at some steady voltage. 

iitemnesults Of the simulation for the three nonlin- 
earities considered are shown in pieure Iii 7: 
Pee oc aling 

When simulating problems on an analog computer, scaling 
is quite important. There are many approaches to this 
feemaaqie, What is of interest here 1s to point out that 
it is necessary to perform magnitude and frequency scal- 
ing on the analog simulation of the nonlinearities simula- 
Eamets study. this 1s true for the two nonlinearities 


with hysteresis. The simulation for the saturation nonlin- 


cate eonly requires amplitude scaling. 


a 








PRCURUMoDe oa ebacklash characteristic curve... 





FIGURE B-3b: Analog simulation for backlash. 
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Cillation input frequency from 2.9hz to 3.0hz 
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FIGURE C-23: Signal waveform showing change of subharmonic 
mode from 1/4 to 1/5. 
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